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Abstract 

■ Our earlier Faddeev three-body study in the -fT - -deuteron scattering length, Aj^-di ls revisited 

. here in the light of the recent developments in two fronts: (i) the improved chiral unitary approach 

OO; 

lO \ to the theoretical description of the coupled KN related channels at low energies, and (ii) the new 



and improved measurement from SIDDHARTA Collaboration of the strong interaction energy shift 
and width in the lowest ET~-hydrogen atomic level. Those two, in combination, have allowed us to 
produced a reliable two-body input to the three-body calculation. All available low-energy K~p 

x 

5^ " observables are well reproduced and predictions for the KN scattering lengths and amplitudes, 

(-7rS)° invariant-mass spectra, as well as for A K - d are put forward and compared with results 
from other sources. The findings of the present work are expected to be useful in interpreting the 
forthcoming data from CLAS, HADES, LEPS and SIDDHARTA Collaborations. 
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I. INTRODUCTION 



There has been a fair amount of recent interest in the low energy interaction of the K 
with the nucleon, the few nucleon systems, as well as nuclear matter. From the point of view 
of possible quasi-bound states, particularly the KNN ones in nucleus, quite a few model 
discussions have been made during the last decade. Relevant references may be found in 
Refs. [l-23]. The essential ingredient is, of course, in the two-body KN system. The isospin 
= 0, strangeness = — 1 part of this system has been known to couple strongly to the near 
threshold A(1405) resonance. So one naturally expects that this resonance should play an 
important role in characterizing the low energy i^-nuclear system. We would like to refer 
the reader to a recent review by Hyodo and Jido 24j, and references therein on the topics 
stated above from the point of view of the chiral approach. Having said that, we go directly 



to our current goal of studying the K -deuteron scattering length, Apc-d- 

Here we are following a path very similar to the one we took [3|, |4j almost 10 years ago, 



At that time new data on the strong interaction level shift (AExs) and width (Ti s ) of the 
^"-hydrogen atomic level had become available from KEK [25|, |26| that was found, for the 
first time, to be consistent with the low energy K~p scattering data, 

AE* EK = -323 ± 63 (stat) ± 11 (sys) eV ; Tf s EK = 470 ± 208 (stat) ± 100 (sys) eV. (1) 

That measurement allowed resolving the long lasting kaonic hydrogen puzzle, born from the 



previous kaonic hydrogen X-ray measurements |27H29| leading to an attractive shift of the 
kaonic hydrogen atom Is state, impossible to reconcile with available theoretical approaches 
(see e.g. Ref. |3o|). 



On the theoretical side, inspired by an earlier work 11, a realistic and reliable description 



of the low energy coupled KN system became available 



2] , which was based on the non-linear 



chiral Lagrangian for the interaction of 
hereafter we call this type of approach 



;he octet pseudo-scalar mesons and octet baryons; 
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36] as the Chiral Unitary Model (xUM). This 
nice opportunity was exploited by Kamalov et al. 3t| to perform a three-body calculation 
of the _fT~-deuteron scattering length, within the fixed center (or the scattering length) 
approximation (FCA) in the input KN amplitudes. In the presence of a two-body resonance 
near threshold (here it is the A(1405)) this approximation had been known to lack accuracy. 
So we made a relativistic three-body calculation of this quantity [l] that automatically took 
into account such effects as nucleon-binding, non-staticity (or target recoil), intermediate 



nucleon-hyperon interactions, etc. Also we demonstrated {jj Q] the importance of retaining 
the deuteron D-state, estimated other possible effects not included in the calculation such as, 
three-particle forces, Coulomb interaction, etc., so that the result might well be considered 
as the most comprehensive model prediction for A^-d available then. Our model prediction 
was A K - d = (—1.80 + -zl .55) fm with an uncertainty of about ±10%. Later, a work by 



Sibirtsev et al. 



38 



came out in which an ingenious method was devised in extracting the 
K°d scattering length from the reaction product of pp -^K°K + d 39]. The result appeared 
to prefer a smaller value in magnitude both for the real and imaginary parts: it looked 
as if our model result were off by two to three standard deviations. However, in view of 



some assumptions and approximations adopted in Ref. [38j] (such as the scattering length 
approximation to the final state enhancement factor which might well vary rapidly due to 
the A (1405)), we believed that their prediction be of semi-quantitative nature (we note that 
the quoted values from our model prediction in that work were incorrect). 

In 2005, a new measurement of the i^~-hydrogen atomic level shift and width became 
available from the DEAR Collaboration 40] , 



DEAR 



-193 ± 37 (stat) ± 6 (sys) eV ; T{ 



DEAR 



249 ± 111 (stat) ± 30 (sys) eV, (2) 



but statistically that appeared to be mutually exclusive with the earlier KEK result 25] . 
So within the version of the chiral unitary model used in {3, ^j], we performed a preliminary 
fit to the low energy coupled KN channel observables by including this new information. It 
was found that the resultant x 2 value became acceptable only when some scattering data 
were excluded from the fit, notably all the K~p elastic cross sections! So it looked as if 
the story had gone back again to the pre-1997 period, hence we did not push further our 
investigations on A K - d anew upon being abided by the DEAR data. 

In the meantime, a few new and improved xUM calculation of the coupled KN channels 
observables became available 
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46] . Unlike earlier xU M approaches that included only the 



leading order term, all these models took up to the next-to-leading order (NLO) terms in the 
chiral exp ansion, which were used in the driving terms of the scattering equation. Borasoy 



et al. 
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43] confirmed the consequence of our preliminary observation in a qualitative and 



convincing manner: the inclusion of the DEAR data made some postdicted cross sections 
deviate significantly from the data. Hereafter the paper by Bor asoy , Nissler and Weise 42] 
is referred to as BNW. As for the third model approach 44|, |46[, two types of solution 
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amplitudes were found: one, consistent with the DEAR data while the other, B\ , not 
consistent. 

In order to make clear our point of view toward the DEAR data, two related theoretical 
works should be worth mentioning here. The first one is an ambitious work by Meissner et 
al. 47j in which the following steps were taken to constrain the possible values of A K - d : 



(i) Obtain the equations relating various physical channel KN scattering lengths (o^y), 
and the isospin scattering lengths and a\ with a due account of the isospin breaking 
effects. 

(ii) Use the consequence of (i)to& given experimental value of o,k- p to impose a constraint 
on do and d\ from unitarity. 

(Hi) From the result of (ii) find the range of allowed values for a-^ N other than ax- P - 
(iv) Then finally based on the result of (Hi), find the range of allowed values for A^-d by 
making use of a FCA as found in Ref. 37]. 

With the value of ax- P extracted from the DEAR experiment, almost all the existing 
model results (except perhaps in the case of an old three-body result 48]) for Ax-d were 
shown to stay outside the imposed limit. On the other hand, the corresponding limit on 
the range extracted from the KEK data was shown to be far more accommodating. With 
the accuracy of the model prediction on the limit due to FCA taken into consideration, we 



49] 



interpreted the outcome as the KEK data to be more relevant. A more recent work 
used a chirally motivated separable model to study the coupled KN channels and, at the 
same time, able to directly calculate AE± S and in momentum space. The result showed 
somewhat lower K~p elastic cross sections compared to the data, and the atomic level width 
came out to be quite larger in magnitude than the DEAR data, although the latter was used 
to constrain the fit. Still the authors concluded that the low energy scattering and the DEAR 
data might not be mutually incompatible (note that their improved study ||3o| led to the 
same finding). 

From what we learned from our study stated above, we decided to calculate A^-d with 
no direct recourse to the DEAR data. We first did the calculation by using the models "c" 



and "s" amplitudes from Ref. 



43] as input. Then, for reasons to be explained later, we also 



obtained our own coupled channel KN amplitudes within the improved \U M, and used 
them to obtain A^-d- These will be discussed later. 

Finally, in mid-2011 the SIDDHARTA Collaboration released 51] a more precise mea- 
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surement of AEu and Ti. s , 
ae sidd = _ 2g3 ± 36 ( sta y ± g eV . pS/BD = 541 ± g9 ^ stat ^ ± 22 e y ^ 

To shed a light on the situation with respect to the three sets of data, we summarize below 
the central value and associated total uncertainties 5 tot = 



u stat * u sys 



AE KEK 
AE DEAR 
AE SIDD 



-323 ± 64 eV ; T? S EK = 470 ± 231 eV, 
-193 ± 37 eV ; rf/ Afl = 249 ± 115 eV, 
-283 ± 36 eV ; rf/ DD = 541 ± 92 eV. 



(4) 
(5) 
(6) 



Then, a few comments are in order: 

(%) Comparing DEAR and SIDDHARTA data shows almost identical uncertainties on 
AE\ S and an improvement of about 20% on Ti s , however the discrepancies between the 
central values come out to be about 2.5cr for both AEi s and 

(ii) Discrepancies in the central values between KEK and SIDDHARTA data show agree- 
ments within la for both A£i s and T\ s . 

(in) SIDDHARTA data improve significantly the precision compared to those reported 
by KEK, namely, total uncertainties go down from 20% to 13% for AE\ S and from 49% to 
17% for r ls . 

In consequence, after having included the SIDDHARTA results in the data base, we 
refitted the coupled KN channels amplitudes, then obtained A^-d- With regard to works 
incorporating the SIDDHARTA data, it should be noted that at the two-body level there 
are two new and improved yJJM calculation of the KN amplitudes by Ikeda, Hyodo and 
Weise 16], and by Mai and Meissner 23[. We will discuss those results later. At the 
three-body level, the previous work by Meissner et al. 47J] was improved and came up with 



results 



52j less restrictive than with the input from the DEAR data for A^-d- The inclusion 



of the recoil effect to FCA 



531 ] is expected to modify the prediction somewhat for the better. 



We also mention a recent work 



541 ] on A K - d where a standard separable interaction picture 



as in Ref. 55| was adopted, but with two-body two-channel (KN — ttS) potentials with a pole 
(to mimic the A(1405)) fitted to the principal K~p initiated channels, and the SIDDHARTA 
data. This is a refinement to a somewhat earlier work by the same author 56J. We note yet 
another piece of recent work 57j] in which, as part of input, the SIDDHARTA result was 
incorporated in a FCA calculation of A K - d as well as the corresponding i^~-deuteron p-wave 



scattering volume to predict the p-wave atomic energy shift and width in the i^~-deuterium. 
Two sets of very distinct solutions were reported there. 

The organization of the present paper is as follows: Section [TT] describes the yJJ M ap- 
proach to the coupled KN channels, which is the primary input to our calculation of Ak-<i, 
Section IIII Al embodies a brief description of the three-body calculation we have adopted. 
A detailed account of this subject may be found in Ref. 4|. Numerical method is outlined 
in Section IIIIBI and our results reported and discussed in Section IIII CI Finally, Section IIVI 
serves as drawing our conclusions. We relegate to Appendix |A] some of the useful formulae 
for Section [TlJ 

II. TWO-BODY CHANNELS 

Since the nucleon-nucleon (NN) two-body interaction that we have adopted in the present 
work is the same as that used in Refs. [si, 0], this section is devoted to the coupled KN 
channel interactions that are the primary two-body ingredient to the three-body approach. 
The SU(3) chiral symmetry guided interaction Lagrangian between the pseudo-scalar octet 
meson: <p = (TT,K,r]), and ground state octet baryons: B = N, Y, H (with Y = A, S) 
is what we need to use here. Due particularly to the dominance of the A(1405) resonance 
close to and below the K~p threshold, which is the energy region of our interest, a question 
might naturally arise as to whether one needs to explicitly introduce a bare (or quark core) 
resonance field A(1405) into the description of the meson-baryon Lagrangian? Earlier, such 



an issue was taken up by Veit et al. 



59| within the context of an extended cloudy bag 
model. The conclusion was that the dominant component of A(1405) appeared to be from 
the KN molecular (or bound) state. More recent studies based on the intrinsic five-quark 



components in the A(1405) wave-function 60( or on a large N c expansion argument 24] 



support this scheme. Thus we adopt a chiral interaction Lagrangian which contains only 



the (f) and B fields as introduced above, by simply noting that there is a recent work [57 1 
on the K~al scattering length and scattering volume by using a Lagrangian with explicit 
bare fields such as A(1405). Since no bound state or resonance can be created by any finite 
number of iterations of the input interaction, it is necessary to iterate the interaction to all 
orders, or to use the i^-matrix , N/D methods, etc. In any event, these are equivalent to 
implementing unitarity, hence the name Chiral Unitary Method (xUM) is adopted. 
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A. KN Interactions 



Coupled channel equations determine the reaction among various two-body meson (</>)- 
baryon (B) channels. For our present study these are KN, irY, r]Y, and KE with different 

charge states. By using i, j, k, , etc. as channel indices, the coupled Bethe-Salpeter 

equations for the t-matrices for the scattering process i — > j reads 

Tij = Vij + VikG^Thj, (7) 

k 

where is the transition potential (or the driving term) taken from an effective chiral 
Lagrangian to be discussed below, and Gq is the free meson-baryon propagator for the 
intermediate channel k. We note here that implicit in the above expression are that (i) 
although not essential, the meson-baryon systems are in the two-body center of mass (CM) 
frame, and (ii) the integration is performed over the off-shell four momentum associated 
with channel k. It is also possible to write the above set of equations collectively in a simple 
matrix form 

f = V + VGT, (8) 

where V = {Vij} is the matrix of driving terms, and similarly for T, while G is a diagonal 
matrix with elements Gq in the diagonal. 

Next, in a standard yJJM approach one makes an on-shell ansatz jsj in which both V 
and T are put fully on-shell with respect to the initial CM channel energy, say, W = 
with s being the Mandelstam s-variable for all the channels involved. Then this ansatz 
renders the above equation into an algebraic one for which the solution becomes 

l-VG(s)] 1 V, (9) 

where the diagonal elements of G(s) are now the off-shell momentum integrated Gq(s), often 
called scalar loops and are functions of only s. A small word of caution: our sign convention 



for V and T is the opposite to the one adopted in Refs. 4l|, |42[. So an apparent difference 
in Eqs.(7-9) from the corresponding ones in those publications may be understood. A more 
detailed account is given in Appendix IA1 

Now we need to specify the fully on-shell V (or Vij). As has become standard by now, one 
obtains these driving terms from a Lagrangian with non-linear realization of chiral symmetry. 
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FIG. 1: Diagrams for meson-baryon scattering: seagull or the Weinberg- Tomozawa interaction (a), contact 

interactions or the Next-to-Leading Order (b), direct s-channel (c), and crossed u-channel Born terms. 
Dashed and solid lines are for mesons and baryons, respectively. 
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42j whose notation and description we will closely follow in 



We refer the reader to Refs. 
our subsequent description. 

Four terms from two primary sources are included for V, Fig. [Q as follows: 
(i) The leading order contact term in the Lagrangian providing the commonly called 
Weinberg- Tomozawa (WT) leading order contact term V^ a \ two Born terms: direct s- 
channel and crossed w-channel V^ d \ which contain two axial vector coupling constants 



D and F. As in BNW, we have adopted the frequently used central values [6l|, namely, 
D = 0.80 and F = 0.46, subject to constraint by the nucleon axial- vector coupling constant 
in the chiral limit, i.e., D + F = qa = 1.26. 

(ii) The NLO term is of the second chiral order 0(p 2 ), and contains a chiral symmetry 
breaking measure (Bq) arising from the chiral condensate, as well as u, d, and s current 
quark masses. The dependence on those quantities is eventually converted to that on the 
physical meson masses, using the Gell-Mann-Oakes-Renner relation for the Goldstone boson 
masses. In addition, for our current objective, this part of the interaction depends on seven 
low energy constants [l|: bo, &d, bp, d\, di-, d^ and d^, which should be determined from 
fits to relevant physical observables. As for the first three parameters which are determined 
from our present fit to available low energy cross sections, we will call them the renormalized 



b parameters: bo, bu and bp. The reason for this may be found in Refs. 



16 



171 ]. and in our 



later discussion. For convenience, we also introduce different combinations of driving terms, 
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called WT-, o, s-, and w-model as follows: 



V WT = V {a \ V c = V {a) + V {b \ V s = V {a) + V {b) + V {c \ V u = V {a) + V {b) + V {c) + V {d \ (10) 

so, for example, the s-model contains WT, NLO and s-pole terms. 

For the present three-body study with total strangeness of the input meson-baryon system 
S = —1, we need two separate sets of coupled channels among various two-body physical 
particle states: 

K~ p -»■ K~p , K°n , Avr° , S+tt" , £°7r° , S~tt + , A77 , S°r/ , [K + Z~] , [K°E°] (11) 



K n — ^ K n , A 



7T 



7T 



7T 



(12) 



The first set contains ten channels: N K - p = 10, and the second one six: N K - n = 6. Thus 
for the first set, Eq. fill)) . channels indices z, j for Vij run from 1 to 10. We note that the 
square bracketed channels in the two equations above were not included in our 2002/2003 
works {2,4] although most of the later works have retained them conventionally. As explained 
later, in the present study we began with the [Nk- p = 10, Nx- n — 6] scheme, then based on 
the statistical fit to data, eventually went back to the earlier choice of [N K - p = 8, N K - n = 5]. 



Having specified the channels, the concrete form of the driving term may now 



De identi- 



fied; Refs. 



17 



42[ provide all the necessary information; see also Sec. Ill 3 of Ref. 



24j (note 



the difference in the overall sign and in spinor normalization in these references). These 
driving terms may need to be spin averaged and projected onto the s-wave orbital state 
before fed into the scattering equation. In these expressions of the driving terms, the weak 
decay constant / enters everywhere. Recall that due to the on-shell ansatz its value is 
expected to be different either from the one in chiral limit, or the physical ones for pion, 
kaon, or rj decays. For our objective its value is determined by fits to the scattering data 



whereas in recent works 



16 



171 ] /„. was fixed at the physical value and fx, f v were floated 



around their physical values and determined by fits. 

Once the driving terms are specified, the last quantity we need to take care of, before 
performing statistical fits to determine the scattering t-matrix T, is G(s) whose j th diagonal 
matrix element is a scalar loop, 

d A q i 



GHs) 



(p - <iY 



Mf 



(13) 



ie 



ra- 



ze 
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where p 2 = s, and Mj and rrij are the baryon and meson masses in the j th channel, re- 
spectively. This is a divergent integral, made finite by dimensional regularization which 
introduces subtraction constants aj, j = 1, 2, .... These constants are also to be fixed by 
fits to the scattering data. Since a scalar loop is characterized only by masses of the particles 
involved, assuming isospin symmetry to reduce the number of these subtraction constants 
may be quite relevant, so for the number of channels Nk- p = 10 (8) the number of ajs is 
6 (5). For more on the scalar loops, see Appendix [A] Here it should be useful to mention 
a very recent work [23], which is within the xUM, but without the on-shell ansatz. Fur- 
thermore, (i) it corresponds to the c-model mentioned above, but included more NLO terms 



fsee the reason behind this choice in Ref. 



62j), (ii) only the first six channels {N K - p = 6) 



are included in the fitting procedure that are open in the energy range for the adopted data. 

Before proceeding to the determination of T using statistical fits to data, it should be 
useful to discuss a couple of items. 

First, we need to state where the present work is different from our 2002/2003 study js, J]. 

(%) Our earlier work 4j included only the WT interaction Vwt i n the driving term whereas 
in the present work we go further to include the NLO contribution, so the c-model V c , and 
also the s- model V s . In practice we have excluded the u- model V u for reasons to be stated 
soon below. Note that the form of the WT interaction used then and in the present work 
are equivalent except for a small contribution disregarded in the former. 

(ii) Our earlier work introduced a global form factor in the driving terms to tame the 
otherwise divergent integral for loop functions. We feel that the dimensional regularization 
in the present work is more consistent in spirit with respecting chiral symmetry. 

(Hi) Apart from physical masses, the implementation of the SU(3) symmetry breaking 
was in the different values of the meson decay constant in our previous work, but in the 
values of the subtraction constants in the present study. 

The second item is on the interaction models we have chosen in the present study as just 
stated in (i) above. As mentioned in BNW, upon adopting the on-shell ansatz, and upon 
projecting onto the s (/ = 0) orbital angular momentum state, the w-pole contribution 7^ 
develops logarithmic dependence on the Mandelstam variable s. This type of logarithmic 
dependence is generally mild and the resultant interaction of lesser importance relative to 
other contributions, except that the branch cuts from the higher threshold channels extend 
up to the threshold of some light meson-baryon channels. We should stress that those 
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branch cuts are unphysical originating from the on-shell ansatz. Possible cures may be, 
(a) as in BNW, to eliminate the singularities by matching the amplitude appropriately to 
some constant where relevant, (b) no particular modification to eliminate the branch cuts 
by arguing that they might only affect a less dominant amplitude for elastic tt°A below the 
7rE threshold 45j, or (c]io adopt a static approximation to baryons in the w-pole terms as 



done, for example, in 50]. 



Here in our work, we simply do not include the u-pole contributions in the driving term 
by acknowledging their relatively weak contribution as already stated above. Thus we have 
studied just the c- and s- models. Underlining this is the observation that, as may be easy to 
understand from the static limit, the s— and u— pole terms tend to compensate each other. 
Consequently, a properly regulated u-model is expected to provide resultant amplitudes and 
observables which may well come out somewhere between the corresponding quantities from 
the c— and s— models. We should emphasize here that once the coupled two-body KN 
amplitudes are fed into the coupled three-body equations, the loop momentum integrations 
go down below the threshold of those two-body amplitudes, so it is essential that the result 
of the three-body calculation not be distorted by any unwanted singularities, hence our 
choice above! Since only BNW have presented separately the cases with c— and s— models' 
results, we were compelled to construct our own KN channel amplitudes. 

To be complete, in Table [J we list the physical particle masses 
work. 

TABLE I: Particle masses (in MeV). 



used in the present 



K~ 


K+ 


K° 


p n TV TT + 7T° 




E+ 


S° 


A 




493.68 


493.68 


497.65 


938.27 939.57 139.57 139.57 134.98 


1197.45 


1189.37 


1192.64 


1115.68 


547.51 1321.31 1314.83 
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B. Data base and extraction of adjustable parameters 



The available data concern various channels / observables and suffer from several incon- 
sistencies. Here, we report on the data used to extract the adjustable parameters for our 
two-body amplitudes. 

In the energy range of interest in this work (Pk- < 250 MeV/c), total cross section 
measurements have been performed between 60's and 80's, producing 64H73I] some 90 
data points for the following channels: K~p — > K~p, K°n, An , H + n~, X°7r°, S~7r + . 
Given that the complete data set shows internal inconsistencies, reduced sets are used in 
various fitting procedures and the number of retained data points is not identical in reported 
phenomenological investigations. Here, we have adopted the same approach as in Refs. 3], 



where 52 cross section data 



65 



66 



69 



70 



72] were selected. The method used there 



was to fit all 90 cross section data, as well as the accurate data 



66 



67 



75KT7| for K~p 



reaction rates at threshold, i.e., 

K~p 7T+£- 



7 
R c 

Rn 



2.36 ±0.04, 



K~p 7T-S+ 
K~p — > charged particles 
K~p — > all final states 

K~p 7T°A 



0.664 ±0.011, 
0.189 ±0.015. 



(14) 

(15) 
(~\ f\\ 

K~p — >• all neutral states 
Then, cross section data giving the highest x 2 were removed and the reduced data base 

was refitted. Checking the outcome of several combinations of data sets allowed establishing 

a consistent enough data base. In the present work, in addition, the strong interaction level 

shift (AEi s ) and width (r ls ) of the ^"-hydrogen atomic level data from the SIDDHARTA 



Collaboration 



5l| were also included in the data base, without affecting the acceptable 



consistency of the kept total cross section data. Finally, we added to the data base the 
pion-nucleon sigma term, a n N = 35 ± 10 MeV, which brings in a loose constraint (we will 
come back to this issue at the end of this section). In summary, our data base embodies a 
total of 58 data points. 

Our approach for N K - p = 8 contains 13 adjustable parameters, as reported in Table [TT1 
The uncertainties attributed to our results are those generated by MINUIT, used for 
minimizations. Our c— and s— models lead both to Xdof ~l-2. As much as the subtraction 
constants are concerned, taking into account the uncertainties, we observe that a^^, and 
a V Y, are rather poorly determined. 
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TABLE II: Adjustable parameters in the present work. The weak decay constant / is in MeV. The 
rcnormalized NLO parameters (&o, bp, bp) and the low energy constants in the NLO Lagrangian (di, d,2, 
d 3 , di) are in GeV -1 . The subtraction constants a(/i) are given at [i=\ GeV. 



Parameter 




Model c 


Model s 


/ 






116.2±0.7 


122.4±0.9 


bo 






-0.35±0.01 


-0.40±0.02 


b D 






0.01±0.02 


-0.08±0.01 


by 






-0.02±0.02 


-0.04±0.03 


di 






-0.17±0.02 


-0.11±0.03 


do 






0.07±0.01 


0.06±0.01 


d 3 






0.29±0.01 


0.26±0.02 


di 






-0.34±0.01 


-0.28±0.02 


a KN 


(xl(T 


3 ) 


1.60±0.08 


1.89±0.10 




(xlO" 


" 3 ) 


-8.59±2.33 


0.96±3.29 




(xlO' 


" 3 ) 


1.15±0.33 


-0.05±0.34 


a v A 


(xl(T 


" 3 ) 


-2.30±0.26 


-1.86±0.34 




(xl(T 




2.99±3.91 


-6.35±4.90 


Xd.o.f. 




1.22 


1.21 



Here, we discuss some issues related to the pion-nucleon sigma term a^N- While we want 
to perform our fit by using direct experimental data only, we are interested in the extent to 
which the constraint from the pion-nucleon sigma term may affect the resulting data fit. At 
leading order the sigma term reads, 



a 



7tN 



-2ml(2b + b D + b F ) 



(17) 



It is well known that the pion-nucleon er-term, obtained by Gasser et al. 78[, based 
on nN data analysis and taking into account the current algebra result generated by the 
quark masses, gave <r vr Ar=45±8 MeV. More recent investigations lead either to smaller or 
larger values with respect to the central one, 45 MeV. For instance, (i) a pertubative chiral 



constituent quark model 



79| | finds 55 MeV, an analysis of ttN scattering amplitude via chiral 
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perturbation theory 80] reaches 59 ± 7 MeV, or still a dispersion relations approach 81] . 



using the SAID nN phase-shift analysis 



82] gives cr,rjv=81±6 MeV, (ii) on the other hand, 



lattice QCD calculations lead to 39 ± 4 MeV 



chiral constituent quark models predict 31 MeV 



, 38 ± 12 MeV 
or 37 MeV 



84L 45 ± 6 MeV |85j, while 



In recent works 



44 



46 



49 



8jJ. 



50], including the present one, some of the NLO low energy 



constants (notably bo,t>D and bp) are constrained by using the lowest order terms for meson 
and/or baryon mass formula (Gell-Mann-Okubo mass formula), as well as the pion-nucleon 
sigma term ovjv within the context of the chiral perturbation theory (%PT). 

On the other hand, by calculating the corresponding value from several xUMs for the 



coupled KN system 



43] . one finds lower values of a n N ~ 15 — 30 MeV. However, as 



discussed briefly in Ref . [17] , the latter is the result of the fit to scattering data by iterating 
the driving term (with physical hadron masses adopted) which contains those NLO constants 
to infinite orders. The extracted values of these constants are then expected to be different 
from those obtained by fitting baryon masses perturbatively, see e.g. Ref. 17] . So as 
stated earlier, we denote the former set of the corresponding values as bo, bo, bp, and the 
corresponding sigma term as o"^, as in Ref. 

In the present work, performing minimizations with no constraint on the sigma term, 
CT^jy was found to fluctuate around 15 ~ 20 MeV. Trying to enforce the sigma term in the 
region 50 to 85 MeV resulted in clustering near the lowest limit. While letting that term 
vary in the range 25 to 45 MeV, the minimization went smoother, with little effect on the 
X 2 - So with the loose constraint 25 < < 45 MeV, we performed the final minimizations. 
The parameters reported in Table HI] lead to a® N = 28 MeV and 36 MeV for models c and 
s, respectively 

To end this section, we would like to comment about other channels, K~p — > 77 A, KE, 
with thresholds beyond the energy range investigated in this work. 

With N K - P = 10 fit, the last subtraction constant cike tends to get quite large relative to 
others, with sizeable uncertainty. This indicates that the two KS channels are very likely not 
relevant to the low-energy fit, and the extracted a^ss have no actual substance. To test this 
observation from several N^- P = 10 fits, we dropped the last two KS channels in Eq. pip , 
and ran the effective Nk- v = 10 — 2 = 8 channel models to calculate the corresponding 
observables. The result was found just about ~ 5% different from the original N K -. 
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cases. To confirm this, we also used Nk- p = 8 fit results, then added two extra KS channels 
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from a few Nk- p = 10 fit results, to make an effective N = 8 + 2 = 10 model. Again, 
the latter was found to reproduce the observables calculated in N^- p = 8 configuration 
within ~ 5%. Also with respect to the data set used in the present work, the N K ~ p = 8 
minimizations are slightly better than the ones with Nk- p = 10 in Xdof by about 0.1. So 
we felt justified to adhere to the N K - p = 8 model. Therein, we kept K~p — > r]Y channels 
which have lower thresholds compared to KE final states, but because of lack of low-energy 
data, we could only check the smooth rising of the cross section close to threshold, matching 
the lowest energy K~p — > 77A data points reported by the Crystal Ball Collaboration |s8 1 . 
Finally, recent studies at higher energies can be found e.g. in [88|, |89|] for 77A final state and 



in Refs. 



46 



9fj, [9l| for KE final states. 



C. Results and Discussion 



In this section we present our results and compare them with the data, namely, cross 
sections for K~p initiated reactions, threshold strong channels branching ratios and kaonic 
hydrogen atom Is level shift. Then we will report on our predictions for the scattering 
lengths and amplitudes, as well as for the (irT,) invariant-mass spectra. 

1. Total cross sections 

In the low energy range of our current interest, i.e. P l £ b < 250 MeV/c, the following 
strong coupled channels are open: 

K~ v -> K~p , K°n , Att° , S+tt" , £°7r° , S~vr + , (18) 

and strongly influenced by the 1 = A(1405) resonance below the K~p threshold which 
decays almost exclusively to 7rS. Moreover, while the K°n has a slightly higher threshold, 
all the 7rY channels have lower threshold than that for K~p. 

Our c-model results for all open strong channels are reported in Fig. [2] and compared to 
the data. The s-model gives very close values to those of c-model, so they are not depicted. 

As mentioned above, we have retained 52 data points for those channels. The obtained 
X 2 per data point (x|LJ turns out to be 0.98. Our model allows reproducing the total cross 
section data satisfactorily. This is also the case for the BNW model, also for the c-model 
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(close to the s-one). The only discrepancy between the two models concerns thei^n final 
state, due to the respective fitted data bases. 




FIG. 2: (Color online) Total cross sections initiated by K p. The solid (blue) curves are our fits and the 



dashed (green) lines our results for the BNW model 

dashed curves arc almost identical for 7r + I]~ and ir°Y,° . Experimental data arc from Rcfs. 



42 1 ; both sets are given for the c-model. The solid and 

1QU" 



70 



72|. 



2. Threshold strong branching ratios and kaonic atom 

Table IIHI summarizes theoretical and experimental values for K~p reaction rates: 7, R c 
and R n (Eqs. (14-16)), as well as i^~-hydrogen atomic Is level shift and width. We get 
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=0.99 for the branching ratios and Xd p =0-21 for the SIDDHARTA kaonic hydrogen 
atom data. 

Our c— and s— model results are close to each other and they agree with the data within 
less than la. Results from other calculations, briefly presented in previous sections, show 
similar trends. We will come back to those works in the following section. 



TABLE III: The K~p threshold strong branching ratios, kaonic atom Is level shift and width (in eV). See 
text for experimental data References. Results for BNW models were obtained using their parameters in 
our code. 



Authors [Ref.] 


7 


Rc 


Rn 


-A£i s 


Tu 


Present work (c) 


2.36 


0.646 


0.190 


314 


589 


Borasoy et al. [42] (BNW (c)) 


2.36 


0.655 


0.191 


316 


562 


Present work (s) 


2.40 


0.645 


0.189 


304 


591 


Borasoy et al. [42J (BNW (s)) 


2.27 


0.652 


0.192 


350 


535 


Mai-Meissner [23] (c-type model) 


2.44±0.70 


0.643±0.017 


0.268±0.098 


296±52 


600±49 


Ikeda, Hyodo, Weise [17] (u-modcl) 


2.37 


0.66 


0.19 


306 


591 


Shcvchcnko [54] (one-pole) 








313 


597 


Shcvchcnko [54] (two-pole) 








308 


602 


Cieply-Smejkal [92] (NLO) 


2.37 


0.660 


0.191 


310 


607 


Krcjcirik [931 


2.36 


0.637 


0.178 


296 


761 


Experiment 


2.36± 0.04 


0.664±0.011 


0.189±0.015 


283±36 


541±92 



42] (BNW) was published before 



It is worth recalling that the work by Borasoy et al. 
the SIDDHARTA data release. 

The only significant deviation from data concerns the width I\ a reported in a very recent 
work 93|, based on the first order Lagrangian, solving the Lippmann-Schwinger equation 



with separable interaction potential. 



3. Scattering lengths 

In view of the K~d scattering length investigations, besides K~p initiated processes, 
Eq. (ITT]) , we need to determine amplitudes for the reactions having as initial state K~n, 
Eq. (I12p . for which there are no data. The adopted procedure is then as follows: once the 
parameters for the T-matrices for the K~p channels, Eq. (11 ip . are determined by fit to the 
data, they are used to calculate the T-matrices for the K~n channels, Eq. ( fl2|) . by assuming 
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SU (3) symmetry in the coupling strengths. Those two sets of amplitudes are then used in 
the three-body calculation of A^-d- 

Here along with the corresponding quantity in the K~p initiated channels, we present 
the scattering lengths as given in Table IIVI For " Data" , kaonic atom measurements are used 
to extract the scattering lengths. "Data" for KEK and DEAR are from Weise 94J. For 
SIDDHARTA, we used the improved 95| Deser-Trueman formula to relate the measured 
quantities to the complex K~p — > K~p scattering length 



ae 1s + -ir u 



2afi(ln a — 1)qk- 



with a the fine-structure constant and [i the K -proton reduced mass. 



(19) 



TABLE IV: KN scattering lengths (in fm); a p , a° and a ex are calculated at W = M K - + M p and a n at 
W = M K - + M n . See text for "Data" explanation. 



Model 




a p (K p 




a„(Jf n 


-> K~n) 


<{K°n 


— >K°n) a ex {K p 


-+K°n) 


Present work (c) 




-0.72 H 


-i0.90 


0.86- 


-i0.71 


-0.12 


fi0.90 -1.21 -4 


-i0.37 


Borasoy et al. [42], BNW (c) 




—0.74 H 


-i0.86 


0.61 - 


-i0.71 


-0.24 


-M0.96 -1.09-4 


-i 0.34 


Present work (s) 




-0.69 H 


-i0.89 


0.90 - 


-i0.66 


-0.10 


fi0.87 -1.214 


- iO.38 


Borasoy et al. [42], BNW (s) 




-0.85 H 


-iO.86 


0.49- 


-i0.67 


-0.38 


4- i 1.01 -1.11-4 


-i0.36 


Ikeda et al. [17] (tt-modcl) 




-0.70 H 


-i0.89 


0.57- 


-i0.73 








Mai-Meissner [23] (c-type model) 


(- 


-0.68 ± 0.15) H 


- j (0.90 ± 0.13) 












Shevchenko [54] (one-pole) 




-0.76 H 


-i0.89 












Shevchenko [54] (two-pole) 




-0.74 H 


-i0.90 












" Data" : 


















SIDDHARTA 


(- 


-0.66 ± 0.07) H 


-j (0.81 ±0.15) 












KEK 


(- 


-0.78 ±0.18) H 


- i (0.49 ± 0.37) 












DEAR 


(- 


-0.47 ±0.10) H 


-i (0.30 ±0.17) 













Note that the scattering lengths in Table IIVI have been obtained at the K~p threshold 
(except for the elastic K~n process). In fact, these quantities are very sensitive to the value 
of the threshold at which they are calculated, which is then reflected in the values obtained 
for the A K - d scattering length. These aspects have been discussed in a previous paper [4[]. 

The scattering lengths, Table IIVI show several features as summarized in the following: 
(i) both real and imaginary parts of a p agree with data within la for all models, with the 
only exception being BNW model-c produced before the release of SIDDHARTA data; (ii) 
for the three other scattering lengths, models' predictions for the imaginary parts, as well as 
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for lZe(a ex ) are compatible with each other; (Hi) finally, in the case of lZe(a n ) and 7£e(et°), 
our predictions turn out to be significantly larger than results from other findings. The 
sensitivity of scattering lengths to model ingredients, especially to the NLO contributions, 



are discussed in Ref. 



4- Scattering amplitudes 

It is instructive to investigate the real and imaginary parts of the KN scattering ampli- 
tudes below threshold: ^/i « 1.432 GeV for K~p, 1.433 GeV for K~n and 1.437 GeV for 
K°n. 

Figures |3] and H] show the amplitudes (/) for elastic scattering channels, obtained from 
our model c, as well as those from BNW. An interesting feature is that for K~p — > K~p 
both real and imaginary part of the scattering amplitudes are very close, as predicted by the 
two depicted models. This however is not the case for the K~n — > K~n process, where the 
real part from our model turns out to be significantly larger than the one in BNW, which 
used DEAR data. 

A behavior common to both models is that the maximum of Im(/^- p ->^- p ) is located 
at y/s wl.416 GeV, so in between the Afl405)_and v /sF^ wl.432 GeV. 



DEAR data were a 
by Hyodo and Weise 



so used in Refs. 



50 



96 



97J. A chiral SU(3) coupled- channel dynamics 



96[, embodying only the dominant Kp — 7rS channels, lead to (much) 



larger amplitudes, in terms of absolute values. In a separable meson-baryon approach, Cieply 
and Smejkal [|3(3], using a two-pole configuration, reported larger real part, but comparable 
imaginary part, though with the maximum around 1.4 GeV. In a recent work based on a 



23| investigated the two-pole 



coupled-channels Bethe-Salpeter approach, Mai and Meissner 
structure and found rather narrow spectrum. 

For the K~n — > K~n scattering amplitudes there are fewer predictions available. In 



addition to that by BNW, shown in Fig.HJ Cieply and Smejkal 50| have reported amplitudes 
with TZe(fK'n^K'n) roughly 50% larger than our values and Tm(fK-n^rK-n) comparable to 
our results, but with a narrower structure. 



19 




FIG. 3: (Color online) Scattering amplitudes for K p channel within model c with real and imaginary 



parts from the present work (blue solid curves) with N K - p = 8, and BNW 



42[ (green dotted curves) with 



N K - p = 10, obtained using our code. 




FIG. 4: (Color online) Same as Fig.[3]but for K~n; N K - n = 5(6) for our model (BNW model). 
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5. 7rS invariant- mass observable 



The (vrS)° invariant-mass spectra data were not fitted in the present work. So, in Fig. [5] 

going back to early 80's, 



we compare our prediction (full curve) to the 7r~S + data [98] 
measured at CERN in the K~p — >■ X + 7r~7r + 7r~ reaction at = 4.2 GeV/c; the agreement 
is reasonable. 




c. m. Energy (GeV) 



FIG. 5: (7rS)° invariant mass spectra for K p — >• 7r S + , 7r + S , 7r°S° as a function of the total center-of- 
mass energy. 



Given that A(1405) formation, decaying to 7rS, is about 30 MeV below K~p threshold, a 
reliable extrapolation oiKN — > 7rS requires data on the shape and location of the invariant 
mass spectra for 7r~£ + , 7r + £~ and 7r°S° channels, for which our predictions are depicted 
in Fig. |5] in full, dashed and dotted curves, res pectiv ely. Predictions for all three channels 



were also reported within chiral approaches 



99 



100[. It is worth noting that those results 
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endorse our predictions, namely, i) the magnitude of the 7r~X + peak is higher than those of 
the other channels, which turn out to have comparable strengths; ii) the peak of the 7r~£ + 
channel is located at lower energy than those of the two other charge states. 



mil 



A recent measurement of the 7r°S° was performed by the ANKE Collaboration 
for P lab = 3. 65 G eV/c via the process pp — > pK + 7T°T,° with rather large uncertainties, as 



discussed in 
tion 



1021 ] . and successfully reproduced within a dy namical c hiral SU(3) calcula 



1031 ] . Very recently, HADES collaboration released data 104 . Il05| for pp — > n^Y^K+p 



at = 3.5 GeV, and A(1405) — > 7r ± S =F channels were extracted. To that purpose, and in 
order to sub trac t contributions from pp — > S°(1385)i^ + p — > ir ± 'E :¥ K + p, the experimentally 



determined 

HADES c 
beam data 



lOq cross sections r atio a pv ^yp n 385) k+ p I Cpp-+ a (i 4n.M k+ p was used. 



aboration results 



1051 ] lead to mass spectra different from the rather old kaon 



depicted in Fig. [51 but also fro m p reliminary measurements with photon 



beams at JLab 



107 



109j and SPring-8/LEPS llOj . Final results using electromagnetic 



probes will hopefully allow checking if the mass spectra depend on the entrance channel, as 
claimed by the HADES collaboration jl05 |. 



Moreover, experimental project at J-PARC 1 111 ] on in-flight (K , n) reaction on deuteron 



is expected to d eepen our understanding of the (7rS)° invariant-mass spectra, via the K d — >■ 
(7rS)°n process 



112 



113]. 



Finally, in order to get in sights into th e quasibound K p states, via 7r S ma ss spectra, 



K absorption in deuterium 



100 



102. 



114j , as well as in 3 He and 4 He j^, [f], Ill5| have been 



investigated theoretically. For a recent review on the nature of A(1405) see Ref. 24]. 



6. Summary on two-body interactions 

The main novelty of the present work with respect to our previous investigations 4] is 
the chiral unitary approach to describe K N interactions, as presented above. 

In order to facilitate the reading of the paper, we summarize the ingredients for other 
two-body channels needed to move to the K~d studies, as reported in detail in Ref. ^]. 

a) A^V interactions: 

For the deuteron (d) channel, we used a relativistic separable potential. All the inter- 
actions considered are of rank-1, which allows correctly reproducing the static parameters; 
namely, the triplet effective range parameters at and rt, the D-state percentage value Pp, the 
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quadrupole moment Q, and the asymptotic ratio rj = Ad /As, the 3 Si phase shift, and also 
the deuteron monopole charge form factor up to about 6 fm~ l . In this model the D-state 
probability is Po = 6.7%. All details can be found in Ref. 116 ] . 
b) nN — P 33 and hyperon-proton interactions: 

Using the same models as in Ref. j^], we have found that these interactions do not have 
a significant effect on the K~d scattering length. We therefore will not bring these partial 
waves in our discussion. 

In conclusion, the obtained KN model and handling of the NN two-body channels lead 
to a reliable enough elementary interactions description to be used in the three-body calcu- 
lation. 



III. K'd SCATTERING LENGTH 



Pioneer work by Toker, Gal and Eisenberg 



oretical investigations 



Ifl 0,0, 



57 



117], some three decades ago, initiated the- 



118Ml2lj on the K d scattering length, in s pite of 



lack of data, though foreseen in a near future. For recent reviews see Refs. 



122 



123]. 



A. Three-body equations for the K d system 

Here, we summarize the three-body equations formalism for the K~d system, in which 
the two-body input described in the previous Section will enter. Then we give the expression 
for the K~d scattering length. 

As the calculations are performed in particle basis, the equations below are given in this 
basis. 

The three-body formalism, where the two-body operators connect two states embodying 
particles which are different, lead to the following system of coupled equations, written in 
operator form: 

X ab {s) = Z ab {s) + J2 Z ac(s)R cc ,(s) X c , h {s), (20) 

c, c' 

with s the three-body total energy. Here, a, b, c and d are the indices which specify the 
particles involved in the spectator and the interacting pair three-body channels. X ab is the 
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transition amplitude between channels a and b, and Z ab is the corresponding Born term. 
The two-body operator R cc < connects two different two-body states labeled as c and d . 

In the particle basis the particles belonging to the various isospin multiplets do have their 
physical masses. The number of 3-body channels is thus considerably enlarged compared 
to the case of isospin basis. So, when using our c- or s- model, the K N interactions must 
include the 8 channels coupled to K~p, namely: K~p, K°n, 7r~E + , 7r + £~, 7r°£°, 7r°A, 77E , 
77A, as well as the 5 channels coupled to K~n: K~n, 7r~£°, 7r°E _ , tt~A, 7/E - . 

Now, we specify the values taken by the channel indices in Eq. ( 120]) . Taking into account 
the deuteron and the 13 (N K - p + N K - n = 8 + 5) two-body inputs in our approach, one must 
consider the three-body channels in the particle basis as given in Table |V] where in the first 
and fourth lines the spectator particles followed (in parenthesis) by the associated pair are 
specified. 

TABLE V: The three-body channels in the particle basis. The second line specifies the isospin of the 
two-body sub-system and the third one the label. 



Channel K~(pn) 


n(K p) 


n(K°n) 


n(7r^S + ) 


n(7r + S 


) n(7r°£°) n(vr°A) 


Isospin 

















Label d 


V\ 


V2 


Q 1 




«3 04 


Channel n{r]Ti ) 


n(rjA) 




P (7T-£°) 


p(7T°£-) 


p(ir~A) p{r]Y,~) 


Isospin 





1 


1 


1 


1 1 


Label [i\ 




2/3 




«6 


Q 7 fl 3 



Starting from the formal equations ( 120 



invariant amplitudes are derived 



124- 



, th e final relativistic equations for the rotationally 



126) : 



Xr aTc (Pa,Pc,s) = Z^ aTc (p a ,p c ; S ) 

dpbPl 



E 

b,T b ;b',r b , 



2e b 



Zi Tb (p a ,p b ; s)Rt b - c »'(a b )Xf biTc (p b ,p c ; s) , (21) 



where a b is the invariant energy of the pair in channel b expressed in the three-body center 
of mass system, c a = (J a , S a , I a ) specifies the conserved quantum numbers of the pair in 
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channel a, and r a = (c a , l a , S a ) specifies the three-body quantum numbers in this channel. 
Labels c and r refer to the spin-isospin variables in a given channel. For example, assuming 
that channel a is composed with particle i as spectator and the pair (jk), we define the 
following quantities: 

• Sf. spin of particle i, 

• Si (= Sj + Sfc), Li, and Jj (= Lj + S^): spin, orbital angular momentum, and total 
angular momentum, respectively, of pair (jk), 

• Sj (= Si + «7j), Zj, and J (= Zj + Sj): channel spin, orbital angular momentum of i 
and (jfc), and three-body total angular momentum, respectively. 

The Born terms matrix in the considered case is a 14 x 14 matrix. However, only a few 
terms are non-zero, as explained in what follows: 

(i) Z p = ( K~(pn) | Go | n(K~p) ): exchange of the p between the deuteron and the (K~p) 
pair, 

(ii) Z n = ( K~(pn) \Go\p(K~n) ): exchange of the n between the deuteron and the (K~n) 
pair, 

(Hi) Zk- = (n(K~p) \G \p(K~n) ): exchange of the K~ between the (K~p) and (K~n) 
pairs, 

(iv) Zj^o = ( n(K°n) \Gq\ n(K°n) ): exchange of theii' between the (K n) pairs 

Of course, we must add the co rresp onding symmetric Born terms, for which the explicit 



expressions can be found in Ref. 



126]. 



In order to obtain the rotationally invariant equations, we only have to antisymmetrize 
the Zj7° = (n(K n) \Gq\ n(K n) ) Born term. This is done as described in Appendix C of 
Ref. [j. 

Concerning the i^iV two-body propagators, they are evaluated as described in Sec. II and 



Appendix A of the present paper, and the d propagator as explained in Ref. 



B. Practical calculation 



To end this Section, we consider the K d scattering length, defined as: 



4]- 



A K - d = - lim 1 X dd , (22) 
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where Xdd is the {J = 1~, I = V — 0) partial amplitude for K~d elastic scattering, evaluated 
at the zero limit for the kaon momentum. 

If we retain the contributions of the d+KN two-body channels, we have a system of 14 
coupled three-body channels (see Table PVT) . After angular momentum reduction, we obtain 
(in the particle basis) a system of 28 coupled equations for J7" = 1~, when including the 
deuteron d, the 8 channels coupled to K~p, and the 5 channels coupled to K~n; see Table 

ED 



The singularities of the kernel are avoided by using the rotated contour method [1271 ]. 
and, after discretization of the integrals, this system is transformed into a system of linear 
equations. In order to solve this system, we use the Pade approximants technique which 
leads to a convergent solution from the successive iterated terms. In practice, we have used a 
diagonal [5/5] Pade (constructed with the 11 first iterates), which was found to be sufficient 
to achieve convergence. Even if the dimension of the matrix to be inverted is rather large, 
this is a sparse matrix because of the limited number of non-zero Born terms, and the Pade 
approximants method is much less time consuming to solve the linear system than the usual 
matrix inversion method. 

TABLE VI: Two-body (L, S, J, T) and three-body (I, S, J) quantum numbers in the particle basis, in the 
case J = 1~, for N K - p = 8 and N K - n = 5. 



channel 


name L 


S 


J 


T 


E 


I 




K~(np) 


d 0,2 


1 


1 





1 


J + 


1 














J - 


1 


n(K-p, K°n, 71--E+, tt+E", 7r°E°, 7r°A, 77E , r;A) 


S/l, 2/2, "1,02,0:3, 04,^1,^2 


1/2 


1/2 





1 


J + 


1 














J - 


1 


p(K~n, 7i--E°, 7r°E- ,7r~A, ^E") 


2/3,05,06,07,^3 


1/2 


1/2 


1 


1 


J + 


1 














J - 


1 



C. Results and Discussion 

The zero limit for the kaon momentum corresponds to W = \/i=1431.95 MeV. As pre- 

n 

viously emphasized |4|, the K N scattering length, and in consequence A K - d , show strong 
dependence on W, due to proximity of the dominant A(1405). It is worth to note that 
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the full three-body results are not subject to excessive sensitivities. Actually, the energy 
dependence in the vicinity of zero limit for the kaon momentum is smeared out due to a 
loop momentum integral. 

Our results for A^-d, an d those from other recent works, are reported in Table IVHl We 
observe that the c- and s-models ioiKN lead to almost identical values for A K - d . 



TABLE VII: K~d scattering length (in fm). 



Authors [Ref.] 


A 


K~ 


d 


Present work (c) 


-1.58 


+ i 


1.37 


Borasoy et al. [42] (BNW (c)) 


-1.59 


+ i 


1.59 


Present work (s) 


-1.57 


+ i 


1.37 


Borasoy et al. [42] (BNW (s)) 


-1.67 


+ i 


1.52 


Doring-Meissner [52] 


-1.46 


+ i 


1.08 


Shevchenko [54J (one-pole) 


-1.48 


+ i 


1.22 


Shevchenko [541 (two-pole) 


-1.51 


+ i 


1.23 


Rcvai [128] (one-pole) 


-1.52 


+ i 


0.98 


Revai [1281 (two-pole) 


-1.60 


+ i 


1.12 


Osct et al. [1291 


-1.54 


+ i 


1.82 


Bahaoui et al. [4] 


-1.80 


+ i 


1.55 



A first comment concerns the fact that our present results, both for real and imaginary 
parts of A K - d , come out about 12% smaller in magnitude with respect to our 2003 values |4j|. 
These changes are due on the one hand to the improved KN model presented in Sec. [Til 
and on the other hand to the since then published SIDDHARTA data 51|. With respect to 
numerical procedure in KN T-matrices calculations, we have carefully checked the conver- 
gence of the integral. In our previous model it was ensured by a form factor in the separable 
form of the interaction, while in the present approach, with no explicit form factor, the 
amplitudes are tempered by dimensional regularization and converge, though more slowly. 

Works by other authors (Table IVIIj) show also the impact of the model used for KN 
interactions and the inclusion or not of the SIDDHARTA data, as discussed below. 

Comparing our results with those of BNW, performed before SIDDHARTA data became 
available, we observe that: (i) Real parts for model c are identical, while for the s-model 
BNW values are in between our present and past {4] results. So, the c-model is less affected 
than the s-model by the SIDDHARTA data, (ii) The imaginary parts, for both c- and 
s-models in BNW are closer to our previous results rather than to those reported in this 
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paper, underlining the signif icant con straints brought in by the SIDDHARTA data. 



All recent results 



52 



54 



128 



129 [reported in Table IVHl were obtained considering the 



SIDDHARTA data, except in Ref. [128 1- 

A new parameterization of the KN — > 7rS potentials was successfully performed by 
Shevchenko 54J and embodied in the author's 56] coupled-channels Faddeev equations in 
Alt-Grassberger-Sandhas form. Results for one- and two-pole structure of A(1405) lead 
to almost identical values (Table IVHj) for A K - d , showing that the scattering length is not 
sensitive enough to the poles structure of the A(1405) resonance. Both real and imaginary 
parts are comparable to our values within less than 10%. Note that in Ref. 54] the NN 
interactions contain no .D-state. 



Revai 128] using a similar three-body formalism, finds slightly larger differences between 
one- and two-pole schemes, especially for the imaginary part of the scattering length, which 
turns out be lower than predictions from all other works quoted in Table IVII1 These features 
might be due to the fact that the two-body channels have been investigated using the 
KEK data rather than the SIDDHARTA results. Actually, Shevchenko has reported such a 
sensitivity considering SIDDHARTA 54] versus KEK 56] measurements. 

Finally, two recent works based on the fixed center app roxim ation (FCA), within a non- 
relativistic effective field theory and Faddeev equations 129], lead to discrepancies in real 
and imaginary parts of about 24% and 44%, respectively. However, a careful study 52] on the 
accuracy of predicted A K ~ d puts forward an "allowed" surface in the [JZe(A K - d ), Tm(A K - d )}- 
plane. Interestingly, all the results shown in Table IVHI fall within a rather com pact sub- 
space of that plane. The forthcoming data from the SIDDHARTA Collaboration 13o| will, 
hopefully, provide sufficient constraints on that physical sub-space. 



IV. SUMMARY AND CONCLUSIONS 

In section[TT]we performed a comprehensive study on the two-body KN interactions at low 
energies via a Chiral Unitary Model including the next to the leading order terms and the 
coupled Bethe-Salpeter equations. The emphasize was put on the recent kaonic hydrogen 



data released by the SIDDHARTA Collaboration 51], which allowed a consistent use of 
low energy data to extract the 13 adjustable parameters in our approach and reproduce 
satisfactorily the fitted data. 
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Using SU(3) symmetry and the scattering amplitudes for the processes 
K~ p — > K~p, K°n, Air , £ + 7r~, S°7r°, S~7r + , A77, £°?7, those for the reactions 
K~n — > K~n, An~ , S°7r~, S~7r°, £"77 were also determined, providing the needed 
complete KN inputs to the K~d system investigations. 

Within the two-body studies, we put forward predictions for entities of interest, namely, 
scattering amplitudes and lengths, as well as the (vrS)° invariant-mass spectra. Results 
from other authors having used the SIDDHARTA data give comparable results e.g. for the 
scattering length a^-p^K-p- Data for (vrS)° invariant-mass spectra has never been used 
to constraint the model s, because of p oor statistics. Forthcomin g data, ho pefully accurate 



enough, using hadronic 



101 



104 



105j or electromagnetic probes 107^ HOj are expected to 



significantly improve our knowledge on the A(1405) — > (ttS)° transitions and hence on the 
KN low energy interactions. This latter is also of paramount importance in pinning down 
whether the antikaon is or not bou nd in few-body nuclear systems, for which advanced 
formalisms has been developed 131 ) . 

The elementary operators obtained in section [IT] were then used in our three-body 
equations approach, embodying the relevant NN two-body channel, namely the deuteron, 
which was described using a relativistic separable potential, with the D-state probability 
P D = 6.7%. 

Our prediction for the K~d scattering length is A K - d = -1.58 + i 1.37. Results from 
various calcula tion s need to be compared to the forthcoming data from the SIDDHARTA 
Collaboration 130 ) . 
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Appendix A: Two-body KN amplitudes 



Since we have closely followed the definition and convention used in BNW [42J, not much 
may be needed to repeat the description of the quantities in the KN equations. So, we will 
simply give some supplementary information. 

(i) Relation to our previous notations 



In case the reader may find it useful, the quantities we adopted in our 2003 work 
denoted as LYS = Lyon-Saclay are related to those in BNW as 



-i] 



Vlys = -47tVbw, (Al) 
T lys = -4nT BNW , (A2) 
Glys — -t—Gbnw, (A3) 

47T 

and the corresponding scattering amplitude reads 

Jlys = — on 2 /-Tlys, (A4) 

where s is the Mandelstam variable of the channel, 
(ii) Two-body meson-baryon scalar loop 

In BNW this function is explicitly given in their Eq. (11). We also used the same sub- 
traction scale value, /x = 1.0 GeV, and m and M are the generic meson and baryon masses, 
respectively. For our practical implementation, we have used a somewhat different form of 
the propagator which is completely equivalent to the form in BNW. We also checked the 
equivalence of our form with the one used in Ref. j^. For convenience, we split G into 
three parts: 

Gbnw(s) = G 1 + G 2 + G 3 , (A5) 

where the constant term is 



Gi = a{ji) + 



1 



327T 2 

with a(/i) the subtraction constant. Then 



(A6) 



and 



G *=^-y- M2 ^> < A7) 



G 3 = -J—AmMVz 2 - lln (-z - Vz 2 - 1), (A8) 
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where 



z = 



s-M 2 - 
2mM 



2 

m 



(A9) 



It might be useful to remind that at the threshold: s = (M + m) 2 , z = 1. Also z = — 1 
at s = (M — m) 2 . Then it can be shown that G is analytic in the whole complex s plane 
with an unitarity cut on the positive real axis starting from the threshold: s = (M + m) 2 . 
An apparent singularity at s = in the above expression, which might derive from the 
kinematics at origin, must be absent. Explicitly we looked at the small \s\ behavior and 
found that 



So, in numerical calculation in the three-body code where the vanishing value of s may 
be encountered, those two contributions need to be calculated together, in order to avoid 
possible spurious divergences. 




(A10) 



(All) 
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